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Abstract 

DFT calculations were carried out on significant Fenugreek seed compounds (1-Methylpyridinium-3-carboxylate MPC, Apigenin-8-C-glucose 

ACG, and 2-lsobutyl-3-methoxypyrazine IBMP) as a green source of ecologically friendly Fe, Al, and Cu metal corrosion inhibitors. Complete 

geometry optimizations were performed by DFT-B3LYP/6-31G* to determine any relationship between the chemical structure and corrosion 

inhibition, mostly on metals' surfaces. Global computational parameters of the inhibitors and thermodynamic Gibbs process of adsorption of metals 

were calculated and used to evaluate each corrosion inhibitor's performance. Our findings showed that MPC has the maximum anti-corrosion 

efficiency across all molecules with a physical adsorption mechanism. It exhibited significant inhibition efficiency with Cu when compared with Fe 

and Al, based on the highest electrophilicity index (w) values compared to other inhibitors and its impact on metals in the following manner: 

Cu>Fe>Al. In regards, the range of inhibitors increased by the following order: IBMP>ACG>MPC. Remarkable corrosion inhibition of MPC is 

demonstrated by its unique high electrophilicity, softness (0), and lowest AEgap on the metal surface. These outcomes are close to the experimental 

work. 
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Introduction 

Metals are integrated into a wide range of industrial products, from steel and iron used in construction materials and car parts to 
aluminum used in containers to titanium used in electrical components. The choice of such metals is linked to their longevity, strength, 
temperature resistance, etc. Examples of certain metals; copper, iron, and aluminum, having their particular characteristics, such as 
copper, which has excellent corrosion resistance, high thermal conductivity. These metals can be prone to corrosion depending on the 
usage environment, including oxygen and certain violent anions like Cl- and SO4? (Amin and Khaled, 2010). This corrosion has a 
detrimental impact on the efficiency of the structures that have been designed on the basis of these metals. Density functional theory 
(DFT) has demonstrated an excellent method in modern physical chemistry due to its capacity to provide specific electron correlation 
results at substantially reduced computational costs (Hohenberg and Kohn, 1964). DFT has offered an excellent foundation for 
designing new standards for rationalizing the prediction, interpretation, and interpretation of many aspects of organic compounds (Wang 
et al., (1999); Udhayakala, (2014)). Sulphathiazole (STZ) has been shown to serve as a potent inhibitor of corrosion of mild steel in 
0.5M HCI, and quantum chemistry calculations indicate that STZ adsorbed mostly as a molecular structure on the metal surface, using 
nitrogen, oxygen and also some carbon atoms as its active sites. Analysis has been conducted by Obot et al. (Obot et al., 2012). The 
group used gravimetric and quantum chemical approaches to test the inhibition performance of STZ onto mild steel in HCl at 303- 
333K. Many corrosion inhibitors such as benzil (4-phenylthlosemicarcazone), Benzoin (BN), benzil (BL), and benzoin (4- 
phenylthiosemicarbazone) (BN4PTSC) are used as corrosion inhibitors in the acidic medium of mild steel. 


* Corresponding author: E-mail: gmazideh @ttu.edu.jo ORCID: https://orcid.org/0000-0003-4494-4284 
Received on January 4, 2021; Accepted on March 16, 2021. 
Jordanian Journal of Engineering and Chemical Industries (JJEC), Vol.4, No.1, 2021, pp. 24-33. 


© The author 


24 


Jordanian Journal of Engineering and Chemical Industries (JJEC! Research Paper Vol.4 No.1, 2021 


Kayadibi et al., used the DFT equations to analyze such corrosion inhibitors' mechanism and effectiveness (Kayadibi et al., 2015). The 
DFT model has emerged as an effective method for researching MPC, ACG, and IBMP compounds as corrosion inhibitors of Cu, Fe, 
and Al metals. Fenugreek (Trigonella Fenugreek-graecum) is one such medicinal plant that has rich phytoconstituents substances like 
alkaloids, steroids, saponins, flavonoids, terpenoids, tannins, and anthocyanin (Yadav et al., (2011)). These components are natural and 
act as a pathogenic attack, single therapeutic agents, and protective agents external effect or combined compounds in drug design 
(Thomas et al., (2006)). The seeds of Fenugreek were found to be effective against Gram-negative bacteria like Escherichia coli and 
Salmonella typhi and Gram-positive bacteria such as Staphylococcus aureus (Al-Mazaideh et al., (2020); Aldal'in et al., (2020)). The 
behavior of the aqueous solution of Fenugreek seeds extract has been studied experimentally as anti-corrosion activity on Aluminum 
(Al) metal in acidic media (1.0M HCI solution) using electrochemical techniques. The extract showed promising results in reducing 
the corrosion rate of Al due to the chemical composing of phytochemical constituents presence as Vitexin and Trigonelline. Three main 
chemical constituents of Fenugreek seeds extract were determined by GC-MS analysis in methanolic extract (Figure 1) (Ladha et al., 
(2015); Karahacane et al., (2017)). 

This study aims to find the systematic method for 
analyzing the potential Fenugreek seed's main 
compounds (1-Methylpyridinium-3-carboxylate 
(MPC), Apigenin-8-C-glucose (ACG), and 2- 


Isobutyl-3-methoxypyrazine (IBMP)) to work as an — o) Í 
inhibitor on the metals' surface interaction in the L / a Í "` 
gas phase using (Density Function Theory-Becke3- ‘coy’ Noon 


parameter Lee-Yang-Parr (DFT-B3LYP)). 
Quantum chemical parameters obtained from this 
calculation and calculated Gibbs function for 
adsorption (AGzas) of the potential inhibitors on Fe, 4-Methylpyridinium-3-carboxylate MPC Apigenin-8-C-glucose ACG —_2-\sobutyl-3-methoxypyrazine IBMP 
Al, and Cu atoms will be used to investigate the 
effect of the studied compounds as green corrosion 
inhibitors for three metals. Achieve comparison of 
theoretical results with experimental data that 
requires developing a common strategy in constructing the inhibitor metal complex models. 


Fig.1 Main chemical constituents of Fenugreek seeds methanolic extract. 


1 Materials and Methods 


Three inhibitors (MPC, ACG, and IBMP) were used to investigate their role as potential corrosion inhibitors for Fe, Al, and Cu atoms 
in the gas phase at 298.15 K. DFT calculations were used in this work using B3LYP functional and 6-31G*(d,p) as basis set by Gaussian 
09 (G09) software (Frisch et al., (2009)). All the inhibitors were optimized by DFT-B3LYP calculations. All the computational works 
were done under no restrictions. The results obtained after optimizing the geometrical structures include quantum parameters of 
inhibitors and the thermodynamic role of adsorption of the metal AG.as. The quantum parameters of the final configuration (Figure 2) 
of the possible inhibitors after a complete optimization and estimation the highest occupied molecular orbital (HOMO) and the lowest 
unoccupied molecular orbital (LUMO) of inhibitors in the gas phase are Eyomo, Etumo, 4E gap, N, 6, X, AN and o (Itawi et al., (2019)). 
The author has done the effects of such calculated parameters on Fenugreek seeds' main compounds to associate the experimental data 
very well with the results (Ladha et al., (2015); Karahacane et al., (2017)). 


1-Methylpyridinium-3-carboxylate MPC 2-Isobutyl-3-methoxypyrazine IBMP Apigenin-8-C-glucose ACG 


Fig. 2 The final configuration of the potential inhibitors after optimization (dark blue (N), gray (C), red (O), white (H)). 
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2 Results and Discussion 
After a geometric optimization of all potential structure inhibitors using the Kohn-Sham method at the DFT level, all quantum chemical 
parameters have been achieved. Table 1 displays Fenugreek extract compounds under investigation as MPC, ACG, and IBMP after 
optimization with their orbitals HOMO and LUMO. 

The FMO (Frontier Molecular Orbital theory) of the inhibitors showed the molecule's adsorbing capacity over the metal surface (Gece 
et al., (2009)). The Exomo indicates a molecule's Table 1 Calculated Eyomo and Exumo of the investigated compounds. 
ability to donate electrons to the appropriate inhibitor Inhibitor Exomo (eV) Emo (eV) AE (eV) 
acceptor with unoccupied metal d-orbitals (Table 1). 


Dipole moment 


H R oR MPC -5.387 -4.071 1.316 15.00 
Conversely, the ELumo indicates its ability to accept 
electrons. Therefore, the lower Eumo value indicates ACG -3.396 1.050 4.446 4.90 
that the molecules have a more remarkable ability to IBMP -8.934 -1.460 7.474 1.99 


accept electrons, leading to a better inhibitor 
efficiency (Li et al., (2009)). Also, by increasing HOMO and decreasing LUMO energy level, the inhibitor's binding capacity to the 
metal surface increases. To determine the type of interaction, the values for Enomo and Eumo for Fe, Al and Cu were compared to those 
calculated for MPC, ACG, and IBMP. Energy gaps (AZ gap) is equal to the differences between Enumo and ELoļmo (Erumo-Enovo) for the 
interplay of metals with inhibitors in given in Table 2. 


Table 2 The calculated E:umo-Enomo gap interaction of metals with the inhibitors. 


Inhibitor Fe* Al* Cu* 
Ecumo:innib-Enomo:re Erumo:re-Enomo:inhib ELumo, inhib - Enomo:al Eumo:ar-Enomo :inhib Ecumo:inhib-Enomo:cu = Etumo:cu—Exomo: inhib 
MPC 3.831 5.236 1.914 5.820 3.655 4.152 
ACG 8.952 3.396 7.035 3.829 8.776 2.161 
IBMP 6.442 8.934 4.525 9.367 6.266 7.699 


* Fe/Eyomo=-7.902eV, Fe/E,umo= -0.151 eV; AI/Enomo=-5.985eV, AI/E,ymo=0.433eV; Cu/Eyomo=-7. 726 eV, Cu/E,ymo=-1.235eV (Lide, (2007)) 


Such computational quantum simulations with an inhibitor design may be a successful approach to estimate the inhibitor's effectiveness, 
regardless of the experimental work results (Karahacane et al., (2017)). Metals will then act as a Lewis acid, whereas the inhibitor MPC 
and IBMP might serve as a Lewis base. Therefore, MPC and 
IBMP use HOMO to activate the reaction of metals with LUMO. 
The Ezumo,Feacu minus Exomo,acc gap quantities, and the 
interaction have a certain amount of ionic character, roughly 
falling between 1.9 to 9eV. It is only possible to predict a 
powerful covalent bond if the ELumo,acc-Exomo Fe,Al,cu difference 
is approximately zero (Lide, (2007); Khalil et al., (2008)). ACG 
functions as a cathodic inhibitor for this case. ACG will serve as 
Lewis acid, while Lewis will be based on metals. Metals use 
HOMO to initiate the response to HOMO of ACG. Due to the 
Ecumo, mpcismp-Exomo,Fe,Al,cu gap values, the association has a 
certain degree of ionic character, roughly falling between 2 to 
9eV. MPC and IBMP serve as anodic inhibitors in this situation. 
The primary responsibility for the chemical reactivity is the 


interaction between the HOMO and the LUMO levels based on Fig. 3 Adsorption of ACG on the metal's surface (ACG—Metal complex). 

the FMO. The energy gap at AEgap is an important parameter and 

shows the inhibitor's responsiveness to adsorption on the metallic surface. Furthermore, the reduction of AE gap is increased. This means 
that the inhibitors’ excitability increases and the potency of the inhibitor are improved (Al-Mazaideh et al., (2018)). The study indicates 
that ACG sees Figure 3, which shows the best potential to coordinate metal vacancy d-orbital links by providing and accepting 
electrons. In the meantime, other universal reactivity parameters have been used to research the efficacy of inhibitors (Arslan et al., 
(2009)). Table 3 shows the global hardness (n), electronegativity (X), global softness (o), a fraction of electron transferred (AN), and 
electrophilicity index (œ) of the inhibitors. 
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In Hard-Soft-Acid-Base (HSAB) term, 
the inhibitors' bonding tendencies 
towards the metal atom can be used. 


Table 3 Calculated global quantum chemical parameters of MPC, ACG and IBMP. 


This term tells that hard acid favors the Parameter MPC ACG IBMP 
coordination of bases with hard 1=-Epomo 5.387 3.396 8.934 
characteristics, and soft acids favor A=-Etumo 4.071 1.050 1.460 
coordinating bases with soft ia eee ses 3:197 
na n=(I-A)/2 0.658 1.173 3.737 
characteristics. Therefore, metal atoms o=1/ 1.519 0.852 0.267 
are classified as weak acids, even AN=(Xreinh) /2MreHinn) 0.077 0.178 0.076 
though soft molecules have a small AN=(Xainn)/2MAarMinh) 0.221 0.124 0.152 
Egap, while complex molecules have a AN=(Xcu-Xinn)/2 (Nou Minn) 0.031 0.255 0.051 
high Egap (Al-Mazaideh et al., (2016). eal cee ee r A 
i yapi ' AIF . 5 i 
The Athree inhibitors’ Egap values ocEu2/2m 15.251 8.555 2685 
(Figure 4) increase in the order: (Xr=4.026, nr=3.875, Xa 3.209, na=2-776, Xa=4480, na=3.250) (Awad er al., (2010) 


IBMP>ACG>MPC. Consequently, for 
metals, inhibitors that are classified as soft bases are the most effective (Li et al., (2009)). Accordingly, compared to many other 
inhibitors, MPC (Fig. 4), which has the least E,,, and the largest softness (o), is expected to have good inhibitor effectiveness. This 
result is verified by calculating the inhibitor's softness (o) and thus its reactivity. Compared with other inhibitors, Table 3 shows that 
MPC has the highest value. Furthermore, the lowest 

hardness (n) of MPC among other inhibitors is AE up(V) 

observed (Table 3). This tendency is the reversal of 

what has been achieved for softness (o). As a result, 

the MPC also with highest softness values (o) and the fe 
lowest hardness value (n) is the best inhibitor. Soft 
molecules are more reactive than complex molecules 
(Al-Mazaideh et al., (2016)). The fraction of 
transferred electrons (AN) is also calculated and ACG 
shows that most of the electronic transition to the 

metals’ surface arises from MPC molecules compared 

to other inhibitors. The electrophilicity index (œ) MPC 


measures the ability of the inhibitor to accept electrons 
from metals. The MPC has the highest electrophilicity 
index (œ), especially on Cu metal compared to the 


l g 
N 
ES 
D 


other inhibitor due to low Erumo. Its metal effect in the m AE gap(eV) 


order: Cu>Fe>Al. This is a good understanding of the Fig, 4 AE gap Of all inhibitors. 
published experiment (Ladha et al., (2015); 
Karahacane et al., (2017)). The computational results 
have shown that MPC adsorbs quickly and forms a thick protective cover to support the metals’ surface. A reaction between nitrogen 
atoms and metal sites helps prevent oxide formation, and the MPC-Metal complex is produced. This observation, therefore, confirms 
that it is highly able to accept metal electrons. Metals, therefore, function as the base of Lewis and Lewis acids as all inhibitors (cathodic 
inhibitors). Inhibitor electrons are accepted by metal as a coordinated bond. Depending on the orientation of the inhibitor's optimized 
structure on the space, an inhibitor can also accept electrons from the metal into back-donating links. As donations and re-donation 
processes strengthen the adsorption of inhibitors to metal's surface, the inhibition efficiency will increase. 
Gibbs free adsorption energy (AGzaas) values of the inhibitors on the metal surface are calculated by using the following equations: 
Cathodic inhibitor: 

AGaas= (inhibitor)-u (M) (1) 
Anodic inhibitor: 


AGaas= (M)-p (inhibitor) (2) 


As (u) references to the chemical potential and it equals to the negative value of electronegativity (X) (Li et al., (2009); Al-Mazaideh 
et al., (2017)), the calculated results of AGzaa; for all inhibitors are given in Table 4. 
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Table 4 Calculated A Gaas values in kJ/mol for the investigated inhibitors. 


Inhibitors Fe Al Cu 
; "E ; sect Anodic Cathodic Anodic Cathodic 
Anodic Inhibitors Cathodic Inhibitors Inhibitors Inhibitors Inhibitors Inhibitors 
MPC 67.829 146.66 24.024 
ACG 173.963 95.134 217.767 
IBMP 112.984 191.813 69.180 


Adsorption of all of the Fe, Al, and Cu inhibitors 
demonstrates that all of them have positive values for 20 | 


AGads aS shown in Figure 5. Most inhibitors act as 
anodic inhibitors of all metals except ACG. The amount 15 
of the value of AGaas depends on whether adsorption is 


m= MPC m ACG IBMP 


chemical or physical adsorption. In physical adsorption, 


10 - 
the value of AGga; is between 0 and -40kJ/mol, while the 
value of AGzas in chemical adsorption is between -80 
and -400kJ/mol (Toliwal et al., (2010)). The suggested a 
= Ie T T 1 


mechanism with all inhibitors is spontaneous chemical k 
=mi ia 
X n o 


adsorption because its value of A Gaas is between -80 and 1 ! 


—40kJ/mol except for MPC with Cu, whereas AGaas is 
within the range of 0 to -40 kJ/mol, and therefore it is 
physical adsorption. This is in line with the 
experimental work recently published (Ladha et al., 
(2015); Karahacane et al., (2017)). 


ANFe ANAI ANCu wfe wWAl wCu 


Fig. 5 Quantum chemical parameter values of all inhibitors. 


Conclusions 


Theoretical calculations provided a good view of the impact of Fenugreek seeds' main compounds (MPC, ACG, and IBMP) as green 
inhibitors, mainly on the corrosion rate of Fe, Al, and Cu metals. Our theoretical calculations exhibited the most excellent adsorption 
of the MPC inhibitor mostly on metal surfaces researched, mainly Cu metal, to make this one of the promising corrosion attacks. MPC 
has shown inhibition efficiency since it has the least ELumo, enabling a better linkage with Cu (more than Fe, Al), behaving as an anodic 
inhibitor. More than the maximum electrophilicity of the MPC compared to the other compounds, the electrophilicity of the molecule 
is measured. All inhibitors show non-spontaneous energy, and MPC has the maximum value of (œ), especially compared to other 
compounds, which recommend that inhibition of corrosive environment is mainly due to MPC. It can be concluded that the MPC is the 
best inhibitor due to its maximum value of softness (o), the highest value of AN, and the highest value of (œ). This comes to fulfilling 
the experimental electrochemical findings of MPC nitrogen atoms' contact to surfaces forming MPC—Metal complex that prevented 
oxide formation. 
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Nomenclature 


ACG =Apigenin-8-C-glucose - 
BL =Benzil - 
BN =Benzoin 


BN4PTSC =Benzoin (4-phenylthiosemicarbazone) 
B3LYP  =Becke, 3-parameter, Lee-Yang-Parr 


AE gap =Energy gap eV] 
DFT =Density Functional Theory -] 
AN =Fraction of electron transferred -] 
Enomo =Energy of highest occupied molecule orbital [eV] 
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Eumo =Energy of lowest un-occupied molecular orbital [eV] 
FMO =Frontier Molecular Orbital theory [-] 
AGaas =Gibbs function for adsorption [kJ/mol] 
G09 =Gaussian 2009 [-] 
HOMO = The highest occupied molecular orbital [-] 
HSAB =Hard-Soft-Acid-Base [-] 
IBMP =2-Isobutyl-3-methoxypyrazine [-] 
LUMO =The lowest unoccupied molecular orbital [-] 

MPC =1-Methylpyridinium-3-carboxylate [-] 

STZ =Sulphathiazole [-] 

o =Global softness [-] 

u = Dipole moment [Debye, D] 
n =Global hardness [-] 

Q =Electrophililcity index [-] 

X =Electronegativity [-] 
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